Multiple-type solutions for multipole interface solitons in thermal nonlinear media 
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We address the existence of multipole interface solitons in one-dimensional thermal nonlinear 
media with a step in the linear refractive index at the sample center. It is found that there exist two 
types of solutions for tripole and quadrupole interface solitons. The two types of interface solitons 
have different profiles, beam widths, mass centers, and stabiUty regions. For a given propagation 
constant, only one type of interface soliton is proved to be stable, while the other type can also 
survive over a long distance. In addition, three types of solutions for fifth-order interface solitons 
are found. 

PACS numbers: 42.65.Tg, 42.65.Jx 



I. INTRODUCTION 

Nonlocality of a nonlinear response is a property ex- 
hibited in many nonlinear optical media. Nonlocal soli- 
tons have been found in nematic liquid crystals [H-Q and 
lead glasses [l-Q theoretically and experimentally. They 
present some novel properties, for instance, the large 
phase shift [l^], self-induced fractional Fourier trans- 
form [ll|, attraction between two dark solitons etc. 
Recently, various types of nonlocal surface solitons [isl - 
[l9l |. for example, multipole surface solitons [T5l|. vor- 
tex surface solitons [Ij], and incoherent surface soli- 
tons [ll] , have been found at the interface between a non- 
linear medium and a linear medium (l3l - [l7j or between 
two nonlinear media [IS . ilQj . The surface solitons propa- 
gating at an interface formed by a nonlinear medium and 
a linear medium are found to be stable only when their 
peaks are less than three • Surface solitons propagat- 
ing at an interface formed by two nonlinear media are 
also found [3 j but dipole surface solitons in such media 
can exist only when optical lattices exist. 

Nonlocal multipole solitons are studied in nematic liq- 
uid crystals [lOl and lead glass [H, [2l| for both bulk 
solitons and surface solitons. In nonlocal bulk media, 
multipole solitons are symmetric, and they are stable if 
they contain fewer than five peaks [13, [2l| . For surface 
multipole solitons, the profiles are asymmetric because 
of the existence of boundary conditions, and they are 
stable when the number of intensity peaks is less than 
three [l^. By comparison, we can qualitatively consider 
surface solitons as half of their corresponding solitons in 
bulk media 22[. For example, surface fundamental soli- 
tons can be regarded as half of dipole solitons in bulk 
media. However, for both bulk solitons and surface soli- 
tons in nonlocal media, only one type of solution exists 
in any case. 

In Ref. [l^, we address the existence of the funda- 
mental and dipole interface solitons propagating at the 



interface between two thermal nonlinear media with dif- 
ferent linear refractive indices. Fundamental interface 
solitons are found to always be stable, and the stabil- 
ity of dipole interface solitons depends on the difference 
in linear refractive index. The boundary force effect [2^ 
plays an important role in the stability of dipole interface 
solitons. It is found that the mass center of the funda- 
mental and dipole interface solitons moves to the part 
with higher linear refractive index as the index difference 
between two media increases [T9| . 

In this paper, we study multipole interface solitons in 
thermal nonlinear media. It is found that there exist two 
types of tripole and quadrupole interface solitons and 
three types of fifth-order interface solitons, respectively. 
The phenomenon of two (or three) soliton types is not 
found in other surface solitons or in bulk solitons. 



II. MODEL OF INTERFACE SOLITONS 

We consider a (l+l)dimensional thermal sample occu- 
pying the region —L<x<L. The sample is separated 
into two parts by the interface at x = 0. All parameters 
for the two parts are same except the linear refractive 
indices. The propagation of a TE polarized laser beam 
is governed by the dimensionless nonlinear Schrodinger 
equation, 

(i) on the left, i.e., —L < x < 0, 



.dq 

dz 



(ii) on the right, i.e., < a; < L, 



.dq 

dz 



15^ 

2dx^ 



+ nq ~ Uciq = 0, 



dx^ 



(1) 



(2) 
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where x and z stand for the normalized transverse and 
longitudinal coordinates, q is the complex amplitude of 
the optical field, n is the nonlinear refractive index, and 
> is the difference in linear refractive index between 
two media. Two boundaries (x — ±L) and the interface 
(x = 0) are thermally conductive. Boundary conditions 



2 



can be described by q{iL) = and n{±L) = 0, and the 
continuity conditions at the interface are q{—0) — q{+0) 
and n{—0) = n{+0). Thus, both q and n are continuous 
at the interface. 

We search for sohton solutions for Eqs. ([!]) and ([2]) 
numericahy in the form q{x,z) = w{x) exp{ibz), where 
w{x) is a real function and b is the propagation constant. 
The details of the model and results for fundamental and 
dipole interface solitons can be found in Ref. [l9[. To 
elucidate the stability of interface solitons, we search for 
the perturbed solutions for Eqs. ([T]) and ([2]) in the form 
q = {w + u + iv) exp(i6z), where u{x, z) and v{x, z) are 
the real and imaginary parts of the small perturbations. 
The perturbation can grow with a complex rate a upon 
propagation. Substituting the perturbed soliton solution 
into Eqs. ([T]) and one can get the linear eigenvalue 
problem. 



1 (fv 

1 d^u 
2d^ 



+ bv — nv, 



— -7-1^ — bu + nu + wAn, 



(-L < X < 0), (3) 



and 



(TV 



1 (Pv 
2dx^ 



+ bv ~ nv + UdV, 



— — bu + nu — UdU + 7«An, 



(0 < a; < L), 



(4) 

where An = —2 j_^G{x,x')w{x')u{x')dx' is the refrac- 
tive index perturbation, the response function G{x, x') — 
{x + L){x' - L)/{2L) for x < x' , and G{x,x') = {x' + 
L){x - L)/{2L) for a; > a;' 



III. MULTIPOLE INTERFACE SOLITONS 

The results of tripole interface solitons are shown in 
Fig. [T] The most interesting feature of tripole interface 
solitons is that there exist two different types of solu- 
tions for some given values of and b. For example, 
when Ud — 0.4 and & = 5, two solutions are shown in 
Figs.[lja) and[TJb), respectively. For the first type of so- 
lution (named type I in this paper) as shown in Fig. [IJa) , 
two intensity peaks are located in the left part of the sam- 
ple with a higher index, and one resides in the right part 
with a lower index. The right peak is much higher than 
the left two peaks because the peak of nonlinear refrac- 
tive index n is located at the right side of the interface. 
For the type-II solution shown in Fig. [IJb), almost all 
three peaks are located in the left part of the sample, 
and the left peak is the highest. Due to the difference in 
the profile, the beam widths (in root-mean-square defini- 
tion) for the two types of tripoles are different as shown in 
Fig. lUc), and they decrease monotonically with increas- 
ing b. However, their powers, defined as P = jgpda; 



and increasing monotonically with increasing 6, are ap- 
proximately equal (relative difference is smaller than 2%) 
as shown in Fig. [Hd). 
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FIG. 1: (Color online) (a) and (b) Profiles for two types of 
tripole interface solitons at Ud = 0.4 and 6 = 5. Dashed- 
dotted lines stand for the interface, (c) Beam width versus 
propagation constant at = 0.1. (d) Soliton power versus 
propagation constant at = 0.1, where the inset is enlarged 
for small values of 6. (e) Mass center and beam width versus 
the index difference at b = 5. (f) Regions of existence and 
stability of tripoles in the b — Ud plane, where shadow areas 
show the stable regions. Points correspond to those cases in 
(c) and (d), and circles correspond to those cases in (e). Red 
(thin) and green (thick) lines stand for type-I and type-II 
tripoles, respectively, in (a)-(f). 

It is known that the interface soliton in our model will 
reduce to the soliton in bulk media when Ud — 0. As Ud 
increases, the tripole interface soliton becomes asymmet- 
ric and shifts to the left part with a higher index. From 
Fig. [TJ^e), only the type-I tripole exists for a small value 
of the index difference, and its mass center (defined as 

Xg — ^^^x\q\^dxl P) and center peak move into the left 
part as Ud increases. When Ud overtakes a certain value 
at a fixed 6, for example, > 0.16 at 5 = 5, there ex- 
ist two types of tripoles [Fig. mf)]. It is interesting to 
note that the right peak of the type-I tripole always re- 
mains in the right part, and its mass center moves back 
to the interface as nd increases sequentially [Figs. [Tfa) 
and[TI^e)]. For the type-II tripole, its right peak moves 
into the left part of the sample, and its mass center shifts 
toward the left and approaches a fixed value when Ud is 



large enough, for example, Xg 



-6.5 at 6 = 5 as shown 



in Fig. [He). The asymptotic feature is induced by the 
boundary force effect (23|, and it is similar to that of 
fundamental and dipole interface solitons fl9( . The mass 
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center of the type-II tripole is proportional to its beam 
width, whereas, the change of mass center for the type-I 
tripole is irrelative with its beam width [Fig. [He)]. 

The existence regions of tripole interface solitons are 
found numerically as 6 > 61 (type I) or & < 62 (type II) 
for a given index difference, where bi and 62 are critical 
propagation constants and are given in Figs. [IIc),[TJd), 
andlljf). In the overlay region (61 < b < 62), which 
increases as rid increases, shown in Fig. WiJ), two types 
of tripoles can exist simultaneously for a fixed Ud- For a 
given propagation constant, there also exists a region of 
Ud in which the two types of tripoles exist simultaneously. 
If 6 = 5, for example, the type-I tripoles exist in Ud < 
0.51 and the type-II tripoles exist in Ud > 0.16, then the 
overlay region is 0.16 < Ud < 0.51 [see the circle symbols 
in Figs. HKe) andHKf)]. 




FIG. 2: (Color online) Soliton intensity profiles vary (a) with 
b a.t Hd — 0.2 and (b) with Ud at b — 5 for type-I solutions 
of tripole interface solitons. Soliton intensity profiles vary (c) 
with b a,t Ud ~ 0.2 and (d) with rid at 6 = 5 for type-II 
solutions of tripole interface solitons. 

Figure [5] shows the intensity peaks of tripole solitons 
vary with b and nd for both type-I and type-II solutions. 
From Figs. [2][a) and [21(c), we can see that the soliton 
power increases with an increase in the constant propa- 
gation at a fixed Ud, which has been shown in Fig. (ijd). 
For a fixed Ud (for example, Ud = 0.2), there exist crit- 
ical propagation constants for both type-I and type-II 
solutions as shown in Fig. [IJf). In Fig. [2Ka), the crit- 
ical propagation constant (61 — 1) is marked, and the 
type-I solution of the solitons disappears when b < bi. 
For type-II solutions, when 6 > 62 = 10, marked in Fig. 
[DJc), one cannot find this kind of solution anymore. If we 
fix the propagation constant, there exists critical Ud for 
both type-I and type-II solutions as shown in Fig. [IJf), 
which also can be explained by Figs. (Hb) andjU^d). Fig- 
ures HKb) and [D[d) present the change in intensity peaks 
versus Ud- For Ud > 0.5, the type-I solutions disappear 
[Fig. mb)], while one cannot find the type-II solutions 
for Ud < 0.16 [Fig. [DJd)]. It is noted that, for type-I 
solutions, the changes in the intensity peaks are obvious 



[Fig. mb)], while the intensity peaks of type-II solutions 
change little [Fig. M.d)]. 

The stability regions of tripole interface solitons are 
found as 6 > 63 (type I) and 64 < 6 < 62 (type II) 
for a given Ud, where 63 and 64 are shown in Fig. [Ijf). 
The letters I (shadow region with vertical lines) and II 
(shadow region with horizontal lines) in Fig.[T](f) indicate 
the stable zones in the b — Ud plane. It is obvious that the 
two types of tripoles can not be stable simultaneously for 
a given propagation constant. In the region 64 < < 62, 
the type-II tripole is stable, but the type-I tripole is not. 

For comparison, tripole solitons in bulk thermal non- 
linear media are stable [1, [2l| , but tripole surface soli- 
tons at the interface between a thermal nonlinear medium 
and a linear medium are unstable [l5|. From Fig. [ijf), 
61 = for type-I solutions when Ud = 0.06 and 62 = 
for type-II solutions when Ud = 0.05. Only the type-I 
tripoles exist for Ud < 0.05, and they are stable almost 
in their whole domain since 63 approaches zero when Ud 
approaches zero. For a very large Ud, it is noted that 
both types of tripole interface solitons have their stabil- 
ity regions, unlike their counterparts in surface solitons. 
For the type-II tripoles, almost all the energy resides in 
the higher- index part, which is similar to the fundamental 
and dipole interface solitons [19] and surface solitons [Tsl - 
[isj . For the type-I tripoles, the right peak stih resides in 
the lower-index medium even for a large Ud [Fig. [HJb)]. 
This unusual feature reasonably cannot be explained only 
by the boundary force effect, since the type-I tripole with 
smaller \xg\ (closer to the center of the sample) should 
get the less equivalent force from the boundaries than the 
type-II tripole. 

The results of quadrupole interface solitons are shown 
in Fig. [HI which are very similar to tripole interface soli- 
tons. There also exist two types of quadrupole interface 
solitons as shown in Figs. [31(a) and[3I^b). The type-I 
quadrupole has three peaks in the left and one peak in 
the right, whereas, the type-II solution has four peaks in 
the left part. The properties of the quadrupoles are the 
same as that of the tripoles, except for the value of the 
parameters. It is worthy to note that there exist stability 
regions, i.e., 6 > 63 (region I with vertical lines) for type- 
I solutions and 64 < 6 < 62 (region II with horizontal 
lines) for type-II solutions, for very large index differ- 
ences. Also, we have observed that there is one peak 
of the type-I quadrupole still residing in the lower-index 
part even for a very large nd [Fig. [2Ja)]. 

To confirm the results of the linear stability analysis, 
we simulate the soliton propagation based on Eqs. ([T]) and 
(|^ with the input condition q{x,z = 0) = w{x)[l + p(x)], 
where w(x) is the profile of the stationary wave and p(x) 
is a random function that stands for the input noise with 
the variance 6'^^^^^ = 0.01. Figure [4] presents propaga- 
tions of tripole and quadrupole interface solitons for both 
type-I [Figs. [la),[lc),[le), and[Ig)] and type-II [Figs. 
[4Kb), [^d), [l^f), and[4Kh)]. As expected, multipole in- 
terface solitons, in their stability regions, survive over 
long propagation distances in the presence of the input 
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tons at different conditions. There exist three types of 
fifth-order interface sohtons. If is small, only one so- 
lution exists with three peaks on the left and two peaks 
on the right (type I) as shown in Fig. [S{a). When 
increases, this type of solution disappears, and other two 
types of solutions exist as shown in Figs. [SKb) andlSJc). 
In Fig. [SJb), there are four peaks on the left and one on 
the right (type II), whereas, in Fig. [SJc), almost all five 
peaks reside on the left (type III). The existence regions 
of the three types of solutions for fifth-order interface 
solitons are shown in Fig. [SKd). The type-I solutions ex- 
ist in the region b > bi [above the red (thin) line], and 
type-II solutions can be found in the region b2 < b < b^ 
[between the two blue (thick) lines] . At the region b < b^ 
[below the green (dashed) line] , one can find the type-Ill 
solutions. It is found that the number of solution types 
is different for different values of Ud and b [Fig. [SJd)]. 
For example, when 6 < 62 or 6 > 63, only one type of 
solution exists. When Ud > 0.36 and 61 < 6 < 64, three 
types of solutions can exist simultaneously. 




FIG. 3: (Color online) (a) and (b) Profiles for two types of 
quadrupole interface solitons at Ud ~ 0.5 and 6 = 3. (c) Beam 
width versus propagation constant at Ud = 0.15. (d) Soliton 
power versus propagation constant at — 0.15, where the 
inset is enlarged for small values of b. (e) Mass center and 
beam width versus the index difference at 6 = 3. (f) Regions 
of existence and stability of quadrupoles in the b — Ud plane. 
The style of each figure is the same as that in Fig[l] 
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FIG. 4: (Color online) Propagations of tripole interface soh- 
tons at 6 = 3 for (a) Ud = 0.05, (b) Ud = 0.15, (e) Ud = 0.1, 
and (f) Ud ~ 0.25. Propagations of quadrupole interface soh- 
tons at fo = 2 for (c) na = 0.05, (d) na = 0.25, (g) na = 0.15, 
and (h) Ud = 0.35. 



noise (the top row). The bottom row in Fig. |4] presents 
propagations of multipole interface solitons in their in- 
stability regions. They experience oscillatory instability 
after propagating over a long distance (> 200 Rayleigh 
distances). This distance is long enough to observe the 
interface solitons in experiments, so one experimentally 
can observe the two types of tripole (quadrupole) inter- 
face solitons at the same Ud and b. 

Figure [S] shows the results of fifth-order interface soli- 



FIG. 5: (Color online) Profiles of fifth-order interface solitons 
at (a) Ud = 0.2, (b) Ud = 0.5, and (c) Ud = 0.5. For all cases 
6 = 3. (d) Regions of the existence of fifth-order interface 
solitons in the b ~ Ud plane. Propagations of fifth-order in- 
terface solitons at (e) na = 0.129, 6 = 0.5 and (f) Ud = 0.5, 
b = 2. 



It is known that the solitons in bulk thermal media are 
unstable when the number of peaks is more than 4 [20l . 
[2l| . However, for the fifth-order interface solitons, there 
exists a stability region, although this region (not given 
here) is very small. One of the major reasons for the 
stable fifth-order interface solitons is the existence of the 
interface. For type-I solutions, one can find the stable 
fifth-order interface solitons when both Ud and b are small 
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[Fig. El^e)]. However, there do not exist stable fifth- 
order interface solitons for type-II and type-Ill solutions. 
Although the two types of solutions are unstable, they 
can propagate a long distance as shown in Fig. [5l^f). In 
addition, for the higher-order (more than five) interface 
solitons, stability regions almost do not exist. 

IV. CONCLUSION 

To conclude, we have studied the properties of mul- 
tipole interface solitons in thermal nonlinear media. It 
is found that there exist two different types of tripole- 
quadrupole interface solitons, and three types of fifth- 
order interface solitons. When the linear index difference 
between two media is small, only the type-I solutions of 
the tripole and quadrupole interface solitons exist, and 
they are stable almost in their whole domain. When 
the index difference is large, two types of tripole and 
quadrupole interface solitons can exist stably in differ- 



ent stability regions, unlike their counterparts in surface 
solitons. It is unusual that the right peak of the type-I 
multipole interface soliton still resides in the lower-index 
medium even for a very large n^. In addition, our model 
can support stable high-order solitons with more than 
four peaks, which is quite different from that in uniform 
media, although there only exists a small stability region 
for the interface solitons with more than four intensity 
peaks. The concept mentioned in this paper can be ex- 
tended to other nonlinear systems. 
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